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During 13 nights of Rayleigh lidar measurements at Urbana, IL in 1984 - 86, thirty- 
six quasi-monochromatic gravity waves were observed in the 35 to 50 km altitude region of 
the stratosphere. The characteristics of the waves are compared with other lidar and radar 
measurements of gravity waves and with theoretical models of wave saturation and 
dissipation phenomena. The measured vertical wavelengths (A^) ranged from 2 to 1 1.5 km 
and the measured vertical phase velocities (c z ) ranged from 10 to 85 cm s-L The vertical 
wavelengths and vertical phase velocities were used to infer observed wave periods (T 0 b) 
which ranged from 100 to 1000 min and horizontal wavelengths (A x ) which ranged from 
70 to 2000 km. Dominant wave activity was found at vertical wavelengths between 2-4 km 
and 7-10 km. No significant seasonal variations were evident in the observed parameters. 
Vertical and horizontal wavelengths showed a clear tendency to increase with T 0 b, which is 
consistent with recent sodium lidar studies of quasi-monochromatic waves near the 
mesopause. An average amplitude growth length of 20.9 km for the rms wind 
perturbations was estimated from the data. Kinetic energy density associated with the 
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waves decreased with height, suggesting that waves in this altitude region were subject to 
dissipation or saturation effects. 


1. INTRODUCTION 

It is now widely recognized that atmospheric gravity waves play a major role in 
determining the large-scale circulation and structure of the middle atmosphere. 
Observational studies of gravity-wave activity have focused on measurements of quasi- 
monochromatic wave parameters and characterization of the continuous gravity-wave 
spectra in order to quantify wave saturation effects and the variability of wave activity with 
altitude, season and location [Fritts, 1984]. Analyses of the temporal frequency and 
vertical wave number spectra of atmospheric wind fluctuations suggest the existence of an 
invariant spectral shape, which is predicted as a consequence of the saturation of vertically 
propagating gravity waves [Dewan and Good, 1986; Smith et al., 1987]. Radar and lidar 
are the predominant remote sensing techniques used to study gravity waves in the middle 
atmosphere. Radars have been successfully used to measure the vertical wave number and 
temporal frequency spectra of gravity waves [e.g.. Smith et al., 1985; Balsley and Carter, 
1982]. Statistical studies of the variability of wave motions and wave characteristics as a 
function of altitude and time have been accomplished in the 60-1 10 km altitude region using 
several MF radar techniques [Meek et al., 1985; Vincent and Fritts, 1987; Ebel et al., 1987; 
Manson and Meek, 1988; Reid and Vincent, 1987]. 

The feasibility of studying atmospheric dynamics with lidar techniques based on 
resonant backscatter from sodium atoms and Rayleigh backscatter from the atmosphere is 
well established. Lidar studies of gravity wave activity generally concentrate on 
observations of quasi-monochromatic events, since wave perturbations usually appear as 
individual events in lidar profiles. The University of Illinois at Urbana-Champaign 
(UIUC) sodium lidar has been used to observe the spatial and temporal variations in the 
sodium layer in the 80-105 km altitude region and to study the propagation of gravity 
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waves near the mesopause [Rowlett et al., 1978; Richter et al., 1981; Shelton et al., 1980; 
Gardner et al., 1986], The parameters of monochromatic gravity waves can be inferred 
from sodium density perturbations [Gardner and Shelton, 1985] and used to infer wave 
saturation and dissipation effects [Gardner and Voelz, 1985]. An extensive characterization 
of monochromatic gravity waves appearing in the sodium layer above Urbana, Illinois was 
published recently by Gardner and Voelz [1987]. 

The altitude region from 25-60 km is generally inaccessible with existing MST radar 
and has not been studied extensively. The lack of observational data in this region provides 
the scientific motivation for Rayleigh lidar studies of atmospheric dynamics. The 
observation altitudes for Rayleigh scatter are not limited to regions containing specific 
atmospheric constituents, as with sodium lidar and other resonance fluorescence lidar 
techniques. Rayleigh lidar measurements are impeded only by the presence of aerosols in 
the atmosphere and by system noise. Rayleigh systems have been successfully used to 
measure atmospheric density and temperature in the 30-80 km altitude region as well as to 
study gravity wave activity [Hauchecome and Chanin, 1980; Chanin and Hauchecome, 
1981; Shibata et al., 1986a]. 

In this paper we investigate the feasibility of using Rayleigh lidars to study gravity- 
wave characteristics in the upper stratosphere. We present the results of 13 nights of 
Rayleigh lidar measurements at Urbana, Illinois (40° 10' N, 88° 10’ W), during which 36 
quasi-monochromatic gravity waves were identified and characterized. The lidar used in 
this study was the UIUC sodium lidar which is a relatively low-performance system when 
compared to state-of-the-art Rayleigh lidars. However, the temporal and spatial resolutions 
are sufficient to observe monochromatic waves and to measure the critical wave 
parameters. In Section 2, the theoretical expressions used to characterize the atmospheric 
density response to monochromatic gravity wave perturbations are derived. The data- 
processing technique used to infer the wave spectrum and horizontal wind perturbations 
from the measured density perturbations is discussed in Section 3. Representative examples 
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of observed wave events are illustrated in Section 4. Seasonal variations of the measured 
gravity-wave parameters and the relationships between the parameters are presented and 
compared to previous radar and lidar observations in Sections 5 and 6. Altitude variations 
of rms wind perturbations are examined in Section 7. Rayleigh lidar performance and data 
processing tradeoffs are discussed in the Appendix. 

2. ATMOSPHERIC DENSITY RESPONSE 
The atmospheric density response to gravity-wave motions is governed by the 
continuity equation. By neglecting diffusion, Gardner and Shelton [1985] have shown that 
the density response can be written in the form 

p(z,t) = e-^poCz-Qz) (1) 

where 

p(z,t) = atmospheric density at altitude z and time t; 

p D (z) = steady-state atmospheric density in the absence of wind perturbations and 
t 

<t>(r,t) = J V • v dx (2) 

-oo 

t 

Gz&O “ j v z dx . (3) 


L= x $ + z z is the position vector where x is the horizontal coordinate and z is the 


vertical coordinate. The wind vector is defined as 
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For monochromatic gravity-wave perturbations, the polarization and dispersion relations 
can be used to show that [Gardner and Voelz, 1987] 

9 Z = yH<|> (5) 

under the condition where X z « 4tcH and co « N, where y is the ratio of specific heats, H 
is the atmospheric scale height, N = [(y-l)g/7H]l/2 is the Brunt-Vaisaia frequency, g is 
the gravitational acceleration, X z is the vertical wavelength of the wave and co is the wave 

frequency. In an isothermal atmosphere the steady-state density decreases exponentially 
with altitude (p Q ~ e* 2 ^) so that, by substituting Eq. (5) into Eq. (1) we obtain for gravity- 
wave perturbations 

p(z,t) = e(Y* 1 ) < f > p 0 (z) . (6) 

For our purposes it is most convenient to work with the natural logarithm of the relative 
density perturbations 


r(z,t) = ln(p/p Q ) = (y- 1 )<p = (— ) 0 

yH 


z • 


(7) 


The gravity- wave polarization and dispersion relations can be used to relate the vertical 
wave number spectra and mean-square values of $ and 0 Z to the power spectrum E x (k z ) 
and mean-square value of the horizontal winds <v x 2 > [Miller et al., 1987] 

E x (k z ) = (yHN)2E <}l (k z ) = N2E 0z (k z ) (8) 


<v x 2> = (yHN)2 «J>2> = n2<0 z 2> . 


(9) 
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As a consequence of Eqs. (7) - (9), the gravity-wave spectrum and mean-square wind 
velocity can be expressed in terms of the Fourier transform and mean- square value of r(z,t) 
[Miller et aL, 1987] 


,yHN, 2 <lR(k z )l 
E x (kz) = (— > T 


( 10 ) 


<v x 2> = (ffll) 2 <r2( z>t )> 

Y-l 

where 


Zc+L/2 

r ik z z 

R(kz) = J r(z,t)e dz . 
z c -L/2 


( 11 ) 


( 12 ) 


L is the altitude range of observations and zc is the altitude at the center of the observation 
interval. 

For the case of low-frequency monochromatic waves, exact solutions for 4> and 0 Z 
were derived by Gardner and Shelton [1985]. 


<J> = In 


1 + 


A e P z 

cos(cot-k*l) 

Y-l 


(13) 


0 Z = yH In 


1 + 


Ae^ z 

cos(cot-k»r) 

Y-l 


where 



= wave amplitude 


( 14 ) 
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P = amplitude growth factor (nr 1 ) 
co = wave frequency (s' 1 ) 
k = k x $ + k z z = wave number vector (nr 1 ) 
k x = horizontal wave number (nr 1 ) 
k z = vertical wave number (nr 1 ). 

The corresponding vertical and horizontal winds generated by an unsaturated gravity wave 
are given by [Hines, 1960] 

v z“-^""“" “ AeP z sin(cot-k*i) (15) 

Y-l * x 


v x = 2— - Ae^ z sin(cot-k»r) (16) 

y-l 


The vertical and horizontal wavelengths are, respectively, X z and X x . The linear density 
perturbation is obtained by substituting Eq. (13) into (7) and noting that the wave amplitude 
is small (<10%) 


r(z,t) = (y-l)ln 


A e ^ z 

1 + cos(cot-k*l) 

Y-l 


= Ae^ z cos(cot-k # i) 


(17) 


Because the density perturbation, r(z,t), and horizontal velocity perturbation, v x , are 
proportional and 90° out of phase, their respective power spectra are proportional, and 
Eqs. (10) and (11) result. 

The kinetic energy per unit mass of a quasi-monochromatic wave is defined as 


KE(k,co)= ^|v x | 2 + j 


2 


(18) 
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where in this case v x and v z are the complex amplitudes of the horizontal and vertical wind 
velocities. For low-frequency waves X z « X x , so that 

KE(k,co) ■ ylv x | 2 = i (^i) 2 (Ae^ z ) 2 . (19) 

L *• y-1 


The presence of quasi-monochromatic waves in the lidar profiles of r(z,t) can be 
determined by examining the spectrum. The spectral signature for a wave is computed 
using Eqs. (10), (12), and (17) [Miller et al., 1987]. 


E x ( k )={-(^) 2 ( Ae pz C ) 2 

yl 


cosh(PL)-l 

2((k-k z )2+p2 ) 


sin2(L(k-k z )/2)' 
((k-k z )2 + p2)_ 


( 20 ) 


Equation (20) is plotted in Figure 1 for X z - 7.5 km, AeP z = 0.01, L = 15 km and several 
values of p. The spectral peak occurs at the spatial frequency X z 'l, and the magnitude in 
m3/s2 is 


E x (k =k 2 ) = |(^~) 2 (AeP z C) 2 


sinh(pL/2) 

L (PU2) 


'l 2 L„_, .[si 

- = jKE(Zc) - 


sinh(pL/2)~ p 
(PL/2) . 


( 21 ) 


where KE(zc) (Eq.(19)) is the wave kinetic energy per unit mass at altitude Zq, the center 
of the observation interval. The vertical wavelength of a monochromatic wave can be 
determined by measuring the spatial frequency of the spectral peak. For most of the data 
reported here, PL<1, so that the square of the bracketed form in Eq. (21) is approximately 
1. Thus the wave amplitude and kinetic energy at altitude z c can be computed from the 
magnitude of the spectral peak using Eq. (21). The vertical phase velocity (c z ) of the 
wave is measured directly from the density perturbation profiles. 
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3. EXPERIMENTAL DATA 

The data used in this study were obtained at Urbana, Illinois, by using the Rayleigh 
scatter photocounts from lidar profiles measured during 1984-86 with the UIUC sodium 
lidar system. The UIUC lidar utilizes a flashlamp-pumped dye laser tuned to the sodium 
resonance line at 589 nm. The operational characteristics of the UIUC lidar system are 
described in Gardner et al. [1986]. The parameters of the lidar arc listed in Table 1. 

The UIUC lidar counts backscattered photons from a single laser shot in discrete time 
intervals that conespond to range bins. The data consist of 100 second photocount profiles 
with a range resolution of 150 m obtained by integrating backscattered photocounts from 
750 laser shots. Absolute density was computed by referencing photocounts near 25 km 
altitude to the atmospheric density specified in the U.S. Standard Atmosphere. 

Photocount profiles can be contaminated by Mie scattering caused by the presence of 
atmospheric aerosols in the propagation path of the laser beam. Rayleigh scattering is 
caused by air molecules and by small aerosols whose the radii satisfy the condition r < 
0.03X. Mie scattering is caused by larger aerosols and particulates such as volcanic ash, 
meteoric dust and water droplets in clouds [Cemy and Sechrist, 1980]. Techniques have 
been developed to detect the presence of aerosols and determine regions of pure Rayleigh 
scattering. These techniques utilize two lasers operating at different frequencies to exploit 
the different wavelength dependence of Mie (X’ 1 ) and Rayleigh (X' 4 ) scattering cross 
sections. Typically, these studies show no aerosols above 25-30 km [Chanin and 
Hauchecome, 1981; Philbrick et al. 1987]. In the UIUC lidar profiles, the lowest usable 
Rayleigh scatter altitude is also limited by photomultiplier (PMT) blanking which is 
employed to prevent receiver overload from strong backscatter returns at low altitudes. For 
the measurements reported here, the PMT gain was switched to maximum at an altitude 
near 30 km. Because the Rayleigh signal decreases rapidly with altitude, the maximum 
usable altitude is dictated by system noise. 
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Rayleigh scatter photocount data can be readily converted to density perturbation 
measurements because the range-scaled signal photocounts are proportional to atmospheric 
density. Atmospheric density is estimated from measured photocount data by range scaling 
the signal photocounts and then multiplying the result by a scaling constant 

p(z,t) = ^z 2 [N(z,t)-Ng(t)] (22) 


where N(z,t) is the photocount at altitude z measured at observation time t and Njj(t) is the 
background photocount per range bin. The density perturbation at altitude z and time t is 
given by (Eq. (7)) 


r(z,t) = ln(p/ Po ) = -1 • (23) 

P 0 (z> 

The second equality holds when the density perturbations are small, which is the case for 
gravity wave effects. The mean atmospheric density Pois estimated by averaging the 

density profiles over the whole observation period 


to+T 0 to+T 0 

p 0 (z) - sp- J p(z,t)dt = i- J ^2 z 2 r N(z,t)-N B ( t )] dt 
lo ^ io ^ 


(24) 


where T 0 is the duration of the observation period (typically 6 to 10 hrs). 

The es timat ed vertical wave number spectrum obtained from the Rayleigh photocount 
data contains a signal component Ex(kz) contaminated by shot noise. While the vertical 
resolution of lidar profiles is fundamentally governed by the width of the receiver range- 
gate, the practical resolution limits are determined by shot noise, which is a function of 
laser power, receiver telescope area, integration time and observation range. The expected 
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vertical wave number spectrum obtained by scaling the spatial power spectrum of the 
density perturbations according to Eq. (10) will have the form (see Appendix) 


£ (^) 2 < I R(kz) I 2 > - E x (k z ) + J- (^) 2 jSl (1 + kS) 2 e lVH 
^ y-1 y-1 X>I T 


where Nj is the total signal photocount within the observation range (z 0 , z Q + L) 
comprising a single lidar profile. The first term on the right side of Eq. (25) is the vertical 
wave number spectrum while the second term is the shot noise component. The linear 
saturation theory predicts that the vertical wave number spectrum of horizontal wind 
perturbations has approximately a k z *3 dependence [Dewan and Good, 1986; Smith et al., 
1987] 


E x (kz) 


The proportionality constant relating the right-hand-side of Eq. (26) to E x (k z ) is on the 
order of unity. In principle, the theory of Dewan and Good [1986] allows exponents other 
than exactly -3. However, much of the existing experimental data indicate that the k z 
exponent is very near -3 and this value will be assumed in the subsequent analysis. 
Equation (26) must be greater than the shot noise floor in order to reliably measure wave 
parameters. Shorter wavelengths can only be observed by decreasing the shot noise 
contamination. The shot noise floor is reduced by increasing Nj, the total signal 
photocount in a profile, or by reducing the altitude range of observations, L. However, 
the use of longer integration periods needed to increase Nt reduces the temporal resolution 
and may average out wave events with shorter periods. Decreasing the altitude range of 
observation will exclude longer wavelength events. Equation (25) is derived in the 
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Appendix and is used to predict the spatial resolution of Rayleigh lidar measurements based 
on the model spectrum. 

The procedure used to process Rayleigh photocount data is similar to the method that 
has been used effectively to process the UIUC sodium lidar data [Rowlett et al., 1978]. 
Photocount data are integrated over adjacent observation times to reduce shot noise levels. 
Density perturbation profiles are then calculated from the integrated photocount profiles. 
The vertical wave number spectrum of each profile is obtained by scaling the spatial power 
spectrum of the density perturbation profile according to Eq. (10). Vertical wavelength 
(X z ) and wave kinetic energy (KE^)) are measured from the vertical wave number 
spectrum. Each density perturbation profile is spatially low-pass filtered in the Fourier 
domain to reduce high frequency shot noise. Vertical phase velocity (cz) is then measured 
by observing the phase progression of the waves in the spatially filtered density 
perturbation profiles. Other characteristics of discrete wave events such as observed period 
are inferred from the measured parameters. 

4. GRAVITY WAVE CASE STUDIES 

To illustrate the analysis technique, we now discuss several wave events observed at 
Urbana, Illinois with the UIUC lidar system. Because the Rayleigh signal levels for the 
UIUC Na lidar are relatively weak compared to state-of-the-art Rayleigh lidars, the 
photocount profiles were integrated for 60 minutes to reduce signal shot noise. To reduce 
computation time the vertical resolution was reduced from 150 m to 300 m by summing the 
counts in each pair of consecutive range bins. The 1-hour profiles were then averaged to 
obtain an estimate of the mean atmospheric density profile for the night. Finally, the 
perturbation about the mean density was then computed for each 60-minute profile. The 
altitude range of interest was usually restricted to 35-50 km. 

The spatial power spectrum of each density perturbation profile was computed by 
calculating the magnitude squared of the discrete Fourier transform. A raised cosine 
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tapering function was used to window the spatial profile in order to control sidelobe effects 
in the power spectrum. The vertical wave number spectrum of the horizontal winds 
associated with the gravity waves was then obtained by scaling this spatial power spectrum 
by [yHN/(Y-l)]2 (Eq. (10)). The spatial power spectra obtained during the evening's 
observations were averaged to reduce the variance of the spectral estimate. The average 
vertical wave number spectrum for 13 August 1984 is shown in Figure 2. Energy density 
E x (k z ) is plotted versus vertical wave number V.j/2% from L'l to the Nyquist frequency 
(2Az)’l where L = 15 km is the observation range and Az = 300 m is the range bin length. 
The dashed line is the estimated shot noise floor. A spectral peak near kjfln = 1.3 x 10*4 
cyc/m (X z = 7.7 km) dominates the power spectrum. The signal-to-noise ratio at this peak 
is almost 7 dB. 

The density perturbation profiles measured on 13 August were spatially filtered with a 
cutoff wave number of 2 x 10*4 cyc/m (1 cyc/5 km). Wave-like structures are clearly 
present in the spatially filtered profiles. In particular, the six consecutive profiles 
corresponding to the observation interval 21:00 to 03:00 (LST) shown in Figure 3 exhibit 
a wave structure with downward phase progression. When this wave was strongest, the 
amplitude inferred from the spatial power spectrum was 2.83%. A vertical phase velocity 
of 0.32 m/s is calculated from the slope of the phase progression. Smaller vertical 
wavelengths were also observed in the Urbana data. Figure 4 shows the filtered profiles 
measured on 29 February - 1 March 1984. The spatial cutoff wavelength is 1.8 km. A 
vertical wavelength near 2.4 km and 0.25 m/s downward phase progression is evident in 
the data. 


5. SEASONAL VARIATION OF GRAVITY WAVE PARAMETERS 
Wave-like structures were observed in almost all density perturbation profiles examined 
for this study. Coherent downward phase progression was identified in two or more 
consecutive profiles at least once during each of the 13 observation nights. Previous 
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Rayleigh lidar observations in the 30-70 km altitude range at Haute Provence indicated that 
wave structures are almost always present in the density profiles for 5 to 60 minute 
integration periods [Chanin and Hauchecome, 1981]. For the Urbana data, a total of 36 
monochromatic waves were observed during the 13 nights. Figure 5 shows the average 
vertical wave number spectrum for the 1984 Urbana data. 

Vertical wavelength (Xz), vertical phase velocity (cz), and wave amplitude (AeP z ) are 
measured directly from the lidar data. For each wave, the observed period T 0 b = A^/c z was 
calculated from the measured values of Xj and c z . The horizontal wavelength 
(X x - T ob /T B , where Tg is the buoyancy period) was inferred using the polarization 

relation for low frequency gravity waves. The measured parameters of each gravity-wave 
event identified during the 13 observation nights are summarized in Table 2. 

The measured vertical wavelengths ranged from 2 to 1 1.5 km. The maximum vertical 
wavelength that can be measured is limited by the altitude range of observations, while the 
minimum value is a function of system noise. Vertical wavelength is plotted versus 
observation date in Figure 6. No distinct seasonal variations are evident. Dominant wave 
activity appears to occur at vertical wavelengths near 2-4 km and 7-10 km. The theoretical 
work of Smith et al. [1987] suggests that the dominant vertical wavelength at stratospheric 
heights is approximately 5 km. The altitude of maximum wind amplitude associated with 
2-4 km vertical wavelength events is predicted to be between 40-60 km [Midgley and 
Leimohn, 1966], Because wave events are most likely to be observed at altitudes where 
the wave amplitude is maximum, our data appear to be consistent with Midgley and 
Leimohn's predictions. 

The values of Xz observed in this study correspond well with previous radar and lidar 
measurements. Rayleigh lidar studies with an XeF laser in the 30-65 km altitude range 
over Fukuoka, Japan, indicate a dominant vertical wavelength near 10 km [Shibata et al., 
1986b]. Rayleigh lidar observations in the same altitude region over Haute Provence 
Observatory in France show dominant vertical wavelengths on the order of 8-15 km 
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[Chanin and Hauchecome, 1981]. Sodium lidar studies of the mesopause region (85-105 
km) above Urbana also show vertical wavelengths in the 2-15 km range [Gardner and 
Voelz, 1987]. Radar measurements of X z in the 60-1 10 km altitude region at Saskatoon 

varied from 0.1-40 km with a median value of 10-15 km [Meek et al., 1985; Manson and 
Meek, 1988]. 

Seasonal variations in vertical phase velocity are shown in Figure 7. The vertical 
phase velocities ranged from 0. 1 1 m/s to 0.85 m/s with a mean value of 0.39 m/s. Shibata 
et al. [1986b] found vertical phase velocities ranging up to 0.33 m/s with minimums near 
zero in January and February. This particular trend is not seen in our data. Chanin and 
Hauchecome [1981] reported that c z is a function of altitude with measured values near 

1 m/s at 50-70 km and less than 0.2 m/s below 50 km altitude. Vertical phase velocities for 
gravity waves in the sodium layer (85-105 km) ranged from 0.36 to 1.75 m/s [Gardner and 
Voelz, 1987]. 

The observed period of each gravity wave event was inferred from measurements of the 
vertical wavelength and vertical phase velocity. Values of observed wave periods ranged 
from 100 to 1000 minutes with a mean value of 240 minutes. The lower value is limited by 
the 60-minute integration time associated with the individual profiles. Due to the presence 
of background wind fields, the observed frequency of the wave is the Doppler-shifted 
frequency rather than the intrinsic frequency. For lidar measurements, the observed 
periods will be biased to values larger than the corresponding intrinsic wave periods 
[Gardner and Voelz, 1987]. The seasonal distribution of T c b is shown in Figure 8. 
Gardner and Voelz did not observe wave events in the sodium layer with periods greater 
than 200 minutes in summer. The wave periods of the few summer events seen in the data 
presented here exceed 200 minutes. However, there does appear to be a slight trend 
towards shorter periods in summer. The radar technique used by Meek et al. [1985] for 
mesospheric observations enabled measurements of background winds and the direction of 
wave propagation. These data were used to compute Doppler shifts of the wave and to 
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estimate the intrinsic period T. Wave periods measured by Meek et al. [1985] varied 
between 5 and 500 minutes while intrinsic periods were as long as 900 minutes. The more 
recent radar data reported by Manson and Meek [1988] suggest that the wave periods are 
shortest in summer which is consistent with lidar data reported here and with the results of 
Gardner and Voelz [1987]. The Manson and Meek data show that in summer and autumn 
the majority of the wave periods are shorter than 40 min while in winter the most frequent 
periods are in the 20 to 60 min range. 

For low-frequency gravity waves, the horizontal wavelength associated with a given 
vertical wavelength and wave period is calculated by applying the dispersion relation. The 
inferred horizontal wavelengths ranged from 40 km to nearly 2000 km. Seasonal 
variations in horizontal wavelengths are shown in Figure 9. Because the observed periods, 
rather than intrinsic periods, were used to calculate the horizontal wavelengths, these 
results may contain substantial errors and should be used with caution. It should also be 
emphasized that in computing the horizontal wavelength, the dispersion relation for an 
isothermal atmosphere was used. In the region between 30 to 50 km, the Brunt-VaisUlh* 
period may increase significantly, thus introducing additional uncertainties in the computed 
horizontal wavelength. Sodium lidar studies did not observe wave events with horizontal 
wavelengths greater than 400 km in summer [Gardner and Voelz, 1987]. However, no 
seasonal trends are evident in the data presented here, and the horizontal wavelengths of the 
few summer events seen in this data set exceed 400 km. The radar data of Meek et al. 
[1985], Reid and Vincent [1987], and Manson and Meek [1988] show horizontal 
wavelengths between about 20 and 300 km. Manson and Meek [1988] observed very few 
long wavelength waves (>400 km) in summer. 

The amplitude of each wave event at the mid-point altitude of the observation range is 
determined by measuring the energy density at the spectral peak and using Eq. (21) to solve 
for Ae^* 0 - The seasonal distribution of wave amplitudes is illustrated in Figure 9. Most of 
the calculated wave amplitudes ranged from 1 to 5%. These values correspond to 



horizontal wind velocities of between 5 and 25 m/s. These values are consistent with radar 
observations at higher altitudes [Reid and Vincent, 1987; Manson and Meek, 1988]. The 
amplitudes of monochromatic waves measured near the mesopause with sodium lidar were 
usually 1 to 5% with maximum values ranging up to 18% [Gardner and Voelz, 1987]. 

Atmospheric scale heights were computed for each observation date from the average 
photocount profile for the evening's observations. The seasonal distribution of 
atmospheric scale heights is shown in Figure 1 1 . Measured scale heights ranged from 6.5 
to 7.6 km with a median value of 7 km. Atmospheric scale heights determined from our 
data appear to be greater in spring and summer months. 

The range of parameters measured in this study is consistent with previous 
measurements of gravity-wave parameters obtained using a wide variety of observation 
techniques [e.g., Gardner and Voelz,1987; Reid and Vincent, 1987; and Manson and 
Meek, 1988]. Dominant vertical wavelengths near 2-4 km and 7-10 km appear evident 
from both the relative occurrence of wave events and Fourier analysis of the data. Shibata 
et al. [1986b] appear to have averaged photocount data over altitude to obtain a 1.5 km 
range bin resolution, and the data presented by Chanin and Hauchecome [1981] were 
smoothed over 3.3 km, precluding observations of short wavelength events. Values of c z 
measured in this study appear larger than values typically observed in the same altitude 
region by the other two Rayleigh lidar studies. No seasonal variations of gravity-wave 
activity are apparent from the wave parameters reported in this study, with the exception of 
a possible trend to shorter wave periods in the summer. More observations in the summer 
months are needed to clarify seasonal activity. Our values of c z in winter are clearly greater 
than zero, in contrast to the near zero values reported by Shibata et al. [1986b]. 

6. RELATIONSHIPS BETWEEN GRAVITY-WAVE PARAMETERS 

The vertical phase velocity is plotted versus vertical wavelength in Figure 12. There is 
considerable scatter in the data presented here. A regression curve of the form c z = CX?? 
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where C is the coefficient and p is the slope has been fitted to the data. Since c z and Xj are 
measured parameters, and consequently, both contain errors, the maximum likelihood 
(ML) regression algorithm described in Voelz and Gardner [1986] was used to estimate the 
power-law fit for the data presented in this figure. This method assumes measurement 
errors in the dependent (Cj) and independent (kj) parameters are statistically independent 

and Gaussian distributed. The ML power-law relation obtained from the data is plotted in 
Figure 12. For comparison, the power-law relation deduced from the sodium lidar data is 
also plotted (dashed line) [Gardner and Voelz, 1987]. The correlation between c z and X z is 
very low (0. 14) in our Rayleigh data. Events with wavelengths near 2-3 km generally had 
0.3 m/s phase velocities when observed with 60-minute integration periods. Shorter 
vertical wavelengths (2-3 km) that were observed near the mesopause with sodium lidar 
were associated with phase velocities near 1 m/s, while the longer wavelength (7-10 km) 
events corresponded to lower phase velocities. The relationship between X z and c z are 
considerably different for the sodium and Rayleigh lidar observations. This may be due to 
the fact that the integration periods were substantially different for the two techniques (i.e., 
10 min for sodium lidar and 60 min for Rayleigh lidar). However, neither Chanin and 
Hauchecome [1981] nor Shibata et al. [1986b] observed large phase velocities in this lower 
altitude region with 20-minute integration periods. 

Both vertical and horizontal wavelengths show a strong tendency to increase with 
increasing periods and are plotted versus the observed period for each wave event in 
Figures 13 and 14. The higher correlation of the data presented in Figure 14 can be 
partially attributed to the manner in which X z is computed [Gardner and Voelz, 1987]. 
Sodium lidar studies of monochromatic wave activity also indicate that X z and X x tend to 
increase with increasing values of Tob [Gardner and Voelz, 1987]. The power-law curves 
for the sodium data are plotted as dashed lines in Figures 13 and 14 for comparison. 
Vertical phase velocity is plotted as a function of the observed wave period in Figure 15. 
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Kinetic energy per unit mass for each monochromatic wave is calculated at the middle 
of the observation range using Eq. (21). The estimated kinetic energy is plotted versus the 
vertical wave number kz/2jt in Figure 16. Much of the scatter evident here can be 
attributed to errors due to shot noise. Although the scatter is significant, the ML power- 
law estimates from the data suggest a k z '2 dependence of kinetic energy. The measured 
distribution can be compared with the vertical wave number dependence of kinetic energy 
predicted by the linear saturation theory [Dewan and Good, 1986]. A monochromatic 
wave is expected to become saturated when its amplitude growth reaches the point of 
convective instability. This implies that at saturation the horizontal perturbation velocity is 
approximately equal to the intrinsic horizontal phase speed of the wave. 


I Vvl 


x 1 sat ~ c ix ~ c x ' v o ~ " 


t b 


(27) 


where vo is the mean background wind velocity in the direction of wave propagation, T is 
the intrinsic period of the wave and Tg is the buoyancy period. The KE distribution 
predicted by saturation theory is then obtained from Eq. (19), 


KEsat 




(28) 


The dashed line in Figure 15 is the predicted distribution of KE sa t- It should be noted that 
although the ML power-law fit to the data suggests a k z _ 2 dependence of kinetic energy, 
because of the large scatter of the data (which results in a very low correlation value of 
0.027) our data should not be taken as conclusive evidence of wave saturation. Indeed, 
most of the data points at vertical wavelengths greater than 5 km fall below the dashed 
KE sa t Ene. By comparison, the distribution of kinetic energies observed in sodium lidar 
studies near the mesopause followed a k z *3 dependence [Gardner and Voelz, 1987]. 
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The kinetic energy of each wave is plotted as a function of the horizontal wave number 
kx/2jt in Figure 17 and observed frequency f D b in Figure 18. The power law fits indicate 
a k x *l and fob**^ dependence. The kinetic energy distributions for the sodium lidar data 
are shown with dashed lines in Figures 17 and 18 [Gardner and Voelz, 1987]. 

The correlation coefficients for the Rayleigh data power law relationships were usually 
less than 0.3. In contrast, sodium lidar measurements showed much higher correlation 
between the gravity wave parameters with values typically ranging between 0.5 and 0.9 
[Gardner and Voelz, 1987]. The significant scatter in our Rayleigh data could be attributed 
to the quality of the data or simply the absence of dominant relationships between the 
gravity- wave parameters in the upper stratosphere. Meek et al. [1985] also computed 
power-law relations for X z and X x versus T and T 0 b- A comparison of wave parameter 
relationships obtained from the Rayleigh lidar data, sodium lidar [Gardner and Voelz, 
1987], and radar studies Meek et al. [1985] is shown in Table 3. 

The differences between the observed relationships deduced from the radar and 
lidar techniques may be caused by genuine geophysical differences between the various 
altitude regions and by the particular measurement biases, resolutions, and accuracies that 
are inherent in each remote sensing technique. It is debatable whether either the radar, 
sodium lidar or Rayleigh lidar techniques adequately sample the wave field. The gravity- 
wave parameters reported by Meek et al. [1985], Manson and Meek [1988] and Reid and 
Vincent [1987] were computed from radar measurements of winds. The Adelaide group 
used a dual beam MF Doppler radar to measure radial winds in the 80-100 km altitude 
region. The two beams sampled scattering volumes that were separated horizontally by 
about 25 km. Wave periods were determined by analyzing the temporal frequency 
spectrum of the measured winds. The horizontal wavelength was determined by 
calculating the phase of the cross spectrum of the winds measured in the two beams and 
comparing the result to a theoretical model. Quasi-monochromatic wave effects were 
assumed to be present in the data when the coherence squared statistic of the cross 
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spectrum was significant (i.e., typically greater than 0.68). The Saskatoon group also used 
an MF radar, 3 spaced receiving antenna systems and a correlation technique to measure 
horizontal drifts of the scattering volumes. These drifts are related to the horizontal wind 
velocities. A cross-spectral analysis technique, similar to that employed by the Adelaide 
group, was then used to identify the gravity-wave events and compute their parameters. 
By employing 3 spaced receiving antenna systems, the background wind field and wave 
propagation direction can also be measured at Saskatoon which enables the intrinsic period 
of the wave to be calculated. 

The MF radar signals result from sporadic partial reflection and/or scattering by thin 
layered structures in the region from about 60 to 120 km altitude. Because of the sporadic 
character of the scattering mechanism, continuous long-term measurements throughout this 
region are rare. The scattering layers rarely occur at night so that observations are restricted 
to daytime with the majority of the measurements obtained within 3 or 4 hours of local 
noon. Depending on whether the echoes come from localized patches of turbulence or by 
"glints" caused by tilted layers, the measured velocity can be related either to the 
background wind or to that associated with wave motions. These uncertainties as well as 
additional ones due to system and observational limitations in measuring the winds may 
contribute to the uncertainties in the gravity-wave parameters deduced from the radar 
measurements. The spectral technique for determining wave periods is well established 
and the measured periods are probably quite accurate. However, the accuracy of the 
measured horizontal wavelengths is not clear. Reid and Vincent [1987] point out that for 
their radar configuration the cross spectrum phase changes slowly for wavelengths greater 
than 200 km. As a consequence, the horizontal wavelength measurements are not reliable 
for the long wavelength waves. The most accurate measurements are obtained for waves 
with horizontal wavelengths comparable to or smaller than the horizontal spacing between 
the scattering volumes. This general observation also applies to measurements made with 
the Saskatoon radar. Neither the Adelaide nor Saskatoon group has quantified how the 
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wavelength measurements are affected by system noise or system resolution. The 
Saskatoon radar has a 3 km vertical resolution while the Adelaide radar has a 2 km vertical 
resolution. The greatest uncertainty in the radar observations lies with the technique for 
identifying quasi-monochromatic wave motions. For both radars, the value of the 
coherence squared statistic is used to determine when the wind velocities within two 
horizontally separated scattering volumes are coherent and indicative of gravity-wave 
motions. When possible, data from different altitudes are compared to verify that the 
winds are compatible with a gravity-wave interpretation. The radar measurements involve 
a great deal of data manipulation and it is not clear how the value of the coherence squared 
statistic is related to the measurement accuracy or whether this statistic is robust enough to 
eliminate spurious wind perturbations from the data. Although additional analyses are 
needed to clarify the accuracies of the radar observations and to determine if the criteria for 
identifying monochromatic waves are appropriate, it is clear that in many cases the radars 
are observing monochromatic gravity waves and the calculated parameters are 
representative of these waves. But, there is also reason to believe that some of the scatter 
in the radar data can be attributed to system limitations and measurement noise and perhaps 
to deficiencies in the wave selection criteria. 

The Rayleigh and sodium lidar techniques also have limitations. Both techniques 
measure density perturbations rather than winds. The sodium and atmospheric density 
perturbations are caused primarily by vertical wind perturbations and the wind divergence 
(see Eqs. (l)-(3». Quasi-monochromatic waves are identified by computing the spatial 
power spectra of the density profiles and then looking for spectral peaks which are 
characteristic of monochromatic gravity waves. For sodium lidar data, the power spectra 
are compared to theoretical models which were developed by assuming that the gravity 
waves are undergoing a uniform Doppler shift throughout the sodium layer or are 
propagating normal to the mean flow. Therefore, the waves which are selected must be 
coherent throughout the 85 to 105 km altitude range of the layer and the vertical wavelength 
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must be relatively constant in this region. This approach eliminates waves which are 
Doppler shifted nonuniformly by the mean wind profile and biases the results towards 
waves propagating normal to the mean flow. In contrast, the radar technique only requires 
coherence of the wind perturbations at two points separated horizontally by about 20 to 30 
km. With lidars, the wave period is determined by measuring the vertical phase velocity or 
by calculating the temporal frequency spectrum of the density perturbations at a fixed 
altitude. 

Gardner and Voelz [1987] presented a detailed analysis of the effects of system 
noise on the accuracies of the sodium lidar measured gravity-wave parameters. Because of 
the excellent height resolution of lidar observations, vertical wavelengths can be measured 
with very high accuracy. Because the techniques for measuring the wave periods are 
similar, the accuracies of the radar and lidar measurements of the observed periods are 
probably comparable. The Saskatoon radar has the additional very important advantage of 
being able to measure the intrinsic periods. The Rayleigh signal levels for the data reported 
here are lower than for the corresponding sodium measurements reported in Gardner and 
Voelz [1987], and so the accuracies of the measured wave parameters are expected to be 
poorer. Because the spectral signature of a monochromatic gravity wave in the Rayleigh 
data is not as complex as that for the sodium data, we believe the Rayleigh technique is less 
selective of the waves. However, the wave perturbations must be coherent throughout the 
30 to 50 km observation region and the vertical wavelength must be relatively uniform in 
this region. As a consequence, there will be some filtering of the wave field by employing 
this selection criteria. Although lidars can operate during daytime, all the reported 
observations of gravity waves were conducted at night when the background noise was 
much less severe. 

In summary, the radar and lidar techniques measure very different but equally 
important characteristics of the monochromatic wave field and they employ radically 
different selection criteria to identify the wave events. The range of parameters 
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representing the true wave field is most certainly more expansive than that measured by the 
Rayleigh and sodium lidars and more restrictive than that measured by the MF radars. 
Both techniques are now making important contributions to the study of wave motions in 
the upper atmosphere and fortunately, the ranges of the various wave parameters measured 
by both techniques do overlap. It is now essential for measurements to be made 
sim ultaneously with both types of systems at the same location so that the strengths and 
limitations of each technique can be completely characterized. 

7. RMS WIND VELOCITY AND KINETIC ENERGY DENSITY 
In addition to studying individual wave events, the Rayleigh lidar data can be used to 
investigate the mean square wind velocity perturbations as a function of altitude. The height 
profile of velocity perturbations provides a quantitative measure of gravity wave activity 
and limi ting processes in the atmosphere. A profile of velocity versus altitude can be 
obtained by scaling the mean-square density perturbations which are calculated at each 
altitude by computing a time average of r 2 (z,t) [Miller et al., 1987] 

*o +T o 

<v x 2 (z)> = <p£^) 2 <r 2 (z)> - (j£2i) 2 _1_ J r 2 (z,t) dt (29) 

Y-l Y-l l ° to 

where T 0 is the duration of observation interval. An altitude profile of rms wind velocities 
for Urbana observations on 8 March 1984 is plotted in Figure 19. To produce this plot 
each 60-minute density perturbation profile was first spatially low-pass filtered using a 
cutoff wavelength of 5 km to reduce shot noise contamination. The mean-square density 
perturbation over the evening was then computed for each range bin and multiplied by the 
scaling factor [yHN/(Y-l)] 2 . The oscillations in the data plotted in Figure 19 result from the 
presence of a slow moving dominant wave . The range of rms velocities is 1 to 4 m/s. 
Wave amplitude clearly grows with altitude, but the large oscillations make it difficult to 
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accurately estimate a growth length. The effects due to these long period oscillations can be 
partially eliminated by averaging the mean-square wind velocity over more observation 
dates. The rms wind profile from 35-50 km for eight of the Urbana data sets is presented 
in Figure 20. Only a portion of the Urbana data can be utilized in this calculation due to 
differences in observation altitudes. The average velocity perturbations increase with 
altitude with values near 1.75 m/s at 35 km to 3.5 m/s at 50 km. Since the data were low- 
pass filtered, only gravity-wave activity corresponding to wavelengths greater than 5 km is 
included. The average amplitude growth length for this profile of rms wind perturbations 
is 20.9 km. The velocity amplitudes of gravity waves in an isothermal atmosphere are 
expected to grow exponentially with height as j n order to conserve kinetic energy. 

The theoretical growth length for unattenuated waves in this altitude region is 2H r 14 km. 

Since the measured growth length for the Rayleigh lidar observations is longer than the 
theoretical growth length, we can infer that waves propagating in this altitude region are, on 
the average, subject to effects that limit the amplitude growth with altitude. Amplitude 
growth lengths ranging from 10-100 km with a mean value of 19 km were measured for 
monochromatic gravity-wave events observed in the sodium layer [Gardner and Voelz, 
1987]. Sodium lidar observations of the rms wind perturbations at Urbana also show an 18 
km amplitude growth length [Senft et al., 1987]. The wind velocities plotted in Figure 19 
can be extrapolated to a value approaching 20 m/s near 85 km by using the measured 
amplitude growth length of 20.9 km. Velocity perturbations of 20 m/s are comparable to 
values obtained from sodium lidar measurements [Gardner and Voelz, 1987; Senft et al., 
1987]. 

Both dissipation and saturation processes will attenuate gravity waves as they propagate 
upward in the atmosphere. Dissipation processes are independent of wave amplitude, 
while saturation depends upon wave amplitude. To illustrate that waves observed in this 
study are propagating upward at attenuated growth rates, we compute the altitude profile of 
kinetic energy density from mean-square velocity perturbations using the relation 



where the atmospheric density p(z) is obtained by referencing the lidar profiles to the 
standard atmosphere at 35 km. The average kinetic energy density profile computed from 
eight Urbana observation nights is shown in Figure 21. Kinetic energy clearly decreases 
with altitude. Decreasing kinetic energy with increasing height indicates that either (1) 
wave energy is being dissipated as waves propagate upward or (2) waves are propagating 
upward through vertical gradients in the mean background wind flow [Balsley and Garello, 
1985; Vincent and Fritts, 1987]. The scale height for the KE profile plotted in Figure 21 is 
18.2 km. Vincent and Fritts [1987] calculated kinetic energy scale heights in the 60-110 
km altitude region with a MF radar and obtained values ranging from 8.6 to 16.2 km. The 
kinetic energy density profile shown in Figure 20 is similar to the KE density profile 
obtained by Balsley and Garello [1985] with the Poker Flat MST radar. The Poker Flat 
radar was used to measure the mean-square horizontal wind perturbations (E-W direction) 
at altitudes below 25 km and above 60 km. The rms winds computed from the MST radar 
KE profile are approximately 8 m/s at 60 km and 1.4 m/s at 25 km. The rms velocities 
estimated from the Rayleigh lidar data in the 35-50 km altitude region presented above are 
within these bounds. The values also agree with those estimated by the model 
computations of Andrews et al. [1987]. 

8. CONCLUSIONS 

In this paper, we have shown that Rayleigh lidar is an effective technique for exploring 
the dynamics of the atmosphere at altitudes between 30-50 km. This is the region which is 
presently inaccessible to radar measurements. The propagation of monochromatic gravity 
waves can be observed by examining altitude profiles of density perturbations and the 
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vertical wave number spectrum of the density perturbations. The data examined in this 
study suggest that substantial gravity-wave activity is present in the upper stratosphere. 
Thirty-six monochromatic wave events were identified and characterized during 13 nights 
of observations. Significant features of the data are the dominant vertical wavelengths 
between 2-4 km and 7-10 km and the tendency for both vertical and horizontal wavelengths 
to increase with an increasing wave period. The kinetic energy density associated with the 
wave activity obtained from the data is consistent with the assumption that the gravity-wave 
amplitudes are being limited either by wave saturation or dissipation effects. The estimated 
magnitude of the rms wind velocity and kinetic energy agree well with values extrapolated 
from MST radar data outside this height region as well as values predicted from model 
computations. It should be mentioned that our analysis has been based on the assumption 
that the lidar-measured density fluctuations are caused by velocity fluctuations associated 
with gravity waves. Our results, therefore, should be checked to see if they are consistent 
with the gravity wave assumption. We have demonstrated, at least qualitatively, that they 
are consistent. 

The quality of the data used in this study is relatively poor in comparison to data which 
can be obtained with wavelength-optimized Rayleigh lidars using state-of-the-art 
technology. The temporal and spatial resolutions of the data presented here are inadequate 
to precisely characterize the entire spectrum of quasi-monochromatic wave behavior. 
Nevertheless, the data analysis techniques developed in this study and the initial 
experimental results clearly demonstrate the potential of using Rayleigh lidar for dynamics 
studies. The results already provide estimates of gravity-wave kinetic energy and wave 
parameters in the upper stratosphere, a region that is inaccessible to radars. 
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Table 1 

UIUC Rayleigh Lidar System Parameters 


4 


Laser 

Laser Wavelength 
Pulse Energy 
Pulse Rate 
Range Resolution 
Receiver Area 
Receiver Efficiency 


Flashlamp Pumped Dye 

589 nm 

50 mJ 

7.5 pps 

150 m 

1.2 m 2 


5% 
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Table 3 

Comparison of Rayleigh Lidar, Sodium Lidar and MF Radar Measurements of the 
Power Law Relationships Between Gravity Wave Parameters 


Rayleigh 
Lidar Data! 

Sodium 
Lidar Data^ 

MF 

Radar Data^ 

Altitude Range: 25-55 km 

80-105 km 

60-1 10 km 

c z =0.091(X z )0.92 

c z = 4.09(^)-0-99 


X z = 0.028(T ob ) 0 - 93 

X z = 0.4 (T ob )0.55 

X z = l72(T)-°- 73 

x x = 0.0i6(To b ) L73 

X x = 0.093(Tob) L52 

^x = 9cr 0 b) 0 - 68 

X x = 39(T)°- 24 

Cz = 31.4(T ob )-0-86 

C Z = 10.2(T ob )-0-54 


KE = 1.3 x lO'^kz)* 1 - 87 

KE = 1.1 x 10* 9 (k z )' 2 - 95 


KE = 4.0x lO^Ckz)* 101 

KE = 5.3 x lO" 4 ^* 1 - 05 


KE = 5.2 x lO-^kz)- 1 - 74 

KE = 9 x 10- 5 (kz)- L59 



^This paper 

^Gardner and Voelz, 1987 
^Meek et al„ 1985 
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Table 2 

Measured Gravity Wave Parameters 


Date 

Observation 

Altitudes 

(km) 

Vertical 

Wavelength 

(km) 

Vertical 

Phase 

Velocity 

(ms* 1 ) 

Observed 

Period 

(min) 

Amplitude 1 

(%) 

15 Jan. 1986 

30-42 

10.5 

0.17 

1029 

1.4 



6 

0.25 

400 

4.0 



4.3 

0.14 

494 

2.68 



4.3 

0.3 

239 

8.25 

25 Jan. 1985 

38-48 

3.8 

0.24 

259 

2.24 



3.3 

0.44 

124 

3.14 



2.2 

0.25 

145 

1.96 

13 Feb. 1984 

30-40 

2.4 

0.4 

100 

2.25 

20 Feb. 1984 

33-48 

10.5 

0.46 

380 

1.79 



2.7 

0.28 

161 

1.74 



2.3 

0.28 

137 

1.74 



2.9 

0.25 

194 

0.89 



2.2 

0.32 

114 

1.18 

29 Feb. 1984 

35-50 

7.2 

0.65 

186 

3.42 



3.4 

0.35 

164 

1.39 



3.4 

0.58 

99 

1.27 



2.4 

0.25 

162 

1.27 

8 Mar. 1984 

35-50 

6.9 

0.23 

500 

2.43 



2.2 

0.27 

139 

2.26 

12 Apr. 1985 

35-50 

3 

0.37 

135 

1.55 

22 Apr. 1986 

30-45 

7.5 

0.48 

259 

2.68 


3.7 

0.58 

106 

5.59 

24 Apr. 1984 

35-50 

6.4 

0.52 

205 

2.47 

35-47 

2.9 

0.37 

129 

3.16 



2.2 

0.19 

190 

4.05 

17 Jul. 1984 

35-50 

9.1 

0.64 

137 

2.62 



2 

0.32 

104 

1.09 

13 Aug. 1984 

35-50 

7.7 

0.32 

417 

2.83 

35-45 

2.4 

0.25 

160 

2.08 

30 Sept. 1984 

35-50 

11.5 

0.85 

225 

2.53 


6.1 

0.38 

268 

4.56 



6.1 

0.61 

165 

5.37 



3 

0.14 

357 

2.19 



2.3 

0.22 

174 

2.47 

12 Nov. 1984 

35-50 

5 

0.65 

128 

2.53 


35-45 

3.5 

0.11 

530 

3.43 


*1 % corresponds to a horizontal wind amplitude of approximately 5 ms' 1 . 
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APPENDIX : RAYLEIGH LIDAR PERFORMANCE EVALUATION 


In this appendix we derive the equations which describe the performance of Rayleigh 
lidar systems. The critical measure of lidar performance is the spatial and temporal 
resolution. Signal shot noise contaminates the data and limits resolution. The data 
processing tradeoffs between integration time, observation range and vertical wavelength 
resolution are described. 

The photocount in a Rayleigh lidar system is proportional to atmospheric density and 
inversely proportional to the square of the range. The expected Rayleigh photocount in the 
absence of wave perturbations can be modeled as 


<N(z)> = 


z o N 0 p(z) _ z qNq j z . Zq )/h 
z 2 P(ZcP z 2 


(Al) 


where N(z) is the measured signal photocount and N 0 is the expected photocount at the 
lower altitude z 0 . The atmospheric density is estimated from the measured photocount by 
subtracting the background noise count, range scaling the result and multiplying by an 
appropriate scaling constant For convenience we will assume the background noise is 

negligible so that 

p(z) = ^jZ^N(z). (A2) 

By substituting Eq. (A2) into Eq. (23) the measured value of r(z) becomes 

1 . 


P</ z ) 



(A3) 
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The discrete Fourier transform of r(z) is given by 


R 



m 


^(mAz) 2 N(mAz) 
p 0 (mAz) 


e ikzfnAz Az 


(A4) 


where Az is the receiver range gate width. The expected power spectrum is obtained by 
averaging the magnitude squared of the discrete Fourier transform. 


ri<m 2 >=XX< 


m 1 


£(mAz) 2 N(mAz) 1 l[ ^(lAz)^(lAz) { 

Pc/mAz) 


>e ik 2 (m-l)Az Az (A5) 


Because N(mAz) is a Poisson process, the photocount fluctuations will introduce shot 
noise into the calculated values of IR(k z )|2. Equation (A5) can be simplified by noting 


J <N„><Ni> 

< N(mAz) N(lAz) > ={ (A6) 

| <N i ^> 2 +<N 1) ^> «=m 


so that 


<|R(k z )| 2 > = 



^(mAz) 2 <N(mAz)> ~ 

e ik z mAz Az 

p 0 (mAz) 



+ 


I 


2 4 

E (mAz) <N(mAz)> , - 
Az J 

Po(mAz) 


(A7) 


The first term in Eq. (A7) is the spectrum of the density perturbations while the second 
term is the shot noise floor. 
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The expression for the shot noise floor can be simplified by approximating the 
summation by an integral, 


S 2 (mAz) 4 <N m > 2 
Az 


Zq+ L 


m 


2 

p 0 (mAz) 


- 1 


2 

z Az 


* 0 No 


> ' z o>/Hfc . Az^d + LjH )2e UH 
Nn Z„ 


(A8) 


We note that the shot noise level is inversely proportional to the photocount at zq and 
grows exponentially with L, the altitude range of the observations. Equation (A8) clearly 
shows that the shot noise floor is smaller for observations at lower altitudes and smaller 
observation ranges. By substituting Eq. (A8) into Eq. (A7) and scaling the result 
according to Eq. (10) we obtain an expression for the measured gravity wave spectrum 


l ( fiV<| RW |’>.E (A9) 

2L y-1 1 1 ^ 2L y- 1 N 0 z o 

It is more convenient to express the shot noise floor in terms of the total photocount 
comprising a profile 


v 4 vt 

Nt= / . < N(mAz) > ” J 


ZjN rt °f e (z -ZoVH N( p 

2 Az 


(A10) 


m 


By solving this equation for N 0 and substituting the result into Eq. (A9), the measured 
gravity wave spectrum becomes 




2L 


Y* 1 


| r«| 


+ —0*~> 2cLIH 
^ Y-1 N t z o 


(All) 
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The theoretical spectrum for saturated gravity waves in the middle atmosphere is 
[Dewan and Good, 1986; Smith et al., 1987], 


Ex(kz)r 


N 2 
2k| ' 


(A 12) 


This expression can be used in Eq. (All) to determine the range of vertical wavelengths 
where the energy density is greater than the noise floor. The result is 


’ 8ji 3 

/7.h 2 y 

LN t 

U-0 


K-a 2 


.L/H 


1/3 




(A 13) 


Equation (A 13) defines the spatial resolution of the Rayleigh lidar. Short wavelengths can 
be observed by increasing the total photocount Nj or by decreasing the altitude range of 

observations L. The photocount return is enhanced by increasing the laser power, 
enlarging the receiver telescope area, or lengthening the measurement integration period. 

The equations derived above can now be used to evaluate the expected performance of 
Rayleigh measurements obtained with the UTUC lidar system. Since L and zq are data 
processing parameters, Nt is the only variable in the above equations that is a function of 

system operating characteristics. The total signal photocount comprising a profile is 
computed using Eq. (A 10) where N 0 is determined from the Rayleigh lidar equation. The 
lidar equation gives the expected number of backscattered photocounts from a laser directed 
into the atmosphere [Cemy and Sechrist, 1980]. The expected Rayleigh photocount per 
pulse from a scattering volume Az thick centered at altitude z is given by 




2q R n A (z)Az 
4ttz 2 


(A 14) 
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where, 

X = laser operating wavelength (m), 

El = laser energy per pulse (j), 
h = Planck’s constant (6.63 x 10*34 j-sec), 
c = speed of light (3 x 10*8 m/sec), 

A r = effective receiver area (m2), 

T a = one-way atmospheric transmittance (0.05), 
or= Rayleigh backscatter cross section (m2), 
n A (z) = atmospheric number density at altitude z (nr 3 ) 

Az = range bin width (m). 

The product of the Rayleigh backscatter cross section and atmospheric density can be 
expressed in terms of atmospheric pressure and temperature [Cemy and Sechrist, 1980] 


/x o «P(z) / 589 .4.0117 

or n A (z) = 3.54 x 10 6 =7-r (- ) 

A,(nm) 


where 


P(z) = atmospheric pressure (mbar) at altitude z 
T(z) = atmospheric temperature (K) at altitude z 


(A 15) 


By substituting Eq. (A15) into (A14) and using the system parameters listed in Table 1 we 
can compute the shot noise floor and vertical wavelength resolution as a function of 
integration time for the UIUC lidar system. 

The performance of the UIUC lidar system was evaluated for observations over the 35- 
50 km altitude range. The predicted shot noise floor is plotted as a function of integration 
time in Figure A.l. The noise floor is clearly reduced with increased integration time. The 
vertical wavelength resolution based upon the model gravity wave spectrum is shown in 
Figure A.2. The expected spatial resolution with 60-minute measurement periods is almost 
2 km. 
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Figure Captions 

Figure 1. Simulated vertical wave number spectrum of a monochromatic gravity wave with 
Xj = 7.5 km, AeP z = 1%, L = 15 km and P' 1 = 14 km, 20 km and 40 km. 

Figure 2. Average vertical wave number spectrum measured during 21: 1 8 to 03:42 LST on 
13 August 1984 using the UIUC lidar at Urbana, IL. The altitude range of 
observations is 35 - 50 km and the vertical resolution is 300 km. The dashed 
line is the estimated shot noise level. 

Figure 3. Spatially filtered density perturbation profiles measured on 13 August 1984 
using the UIUC lidar at Urbana, IL. The low-pass filter cutoff wave number 
was 2 x 10 -4 corresponding to a cutoff wavelength of 5 km. The diagonal 
lines indicate the apparent 0.32 m/s downward phase progression of the 7.7 
km wave. 

Figure 4. Spatially filtered density perturbation profiles measured on 29 February - 1 
March 1984 using the UIUC lidar at Urbana, IL. The low-pass filter cutoff 
wave number was 5.5 x 10 -4 corresponding to a cutoff wavelength of 1.8 
km. The diagonal lines indicate the apparent 0.25 m/s downward phase 
progression of the 2.4 km wave. 
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Figure 5. Average vertical wave number spectrum measured at Urbana, IL on 20 
February, 8 March, 17 July, 13 August, 30 September, and 12 November 
1984 using the UIUC lidar. The altitude range of observations is 35 - 50 km 
with a spatial resolution of 300 m and a temporal resolution of 60 minutes. 
The dashed line is the estimated shot noise floor. 

Figure 6. Seasonal distribution of vertical wavelengths. 

Figure 7. Seasonal distribution of vertical phase velocities. 

Figure 8. Seasonal distribution of observed wave periods. 

Figure 9. Seasonal distribution of horizontal wavelengths. The horizontal wavelengths 
were calculated using the gravity wave dispersion relation by assuming that 
the intrinsic wave period was equal to the observed period T ob . 

Figure 10. Seasonal distribution of wave amplitudes. 

Figure 11. Seasonal distribution of atmospheric scale heights. 

Figure 12. Vertical phase velocity versus vertical wavelength. The solid line is the ML 
power-law fit of the form c z = C(X Z )P where p is the slope and C is the 
coefficient. The power-law fit c z = 4.09(X z )-°" obtained with sodium lidar 
observations is indicated with a dashed line [Gardner and Voelz, 1987]. 
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Figure 13. Vertical wavelength versus observed wave period. The solid line is the ML 
power-law fit of the form X z = C(T ob )P where p is the slope and C is the 
coefficient, the power-law fit X z = 0.4(T ob )*9-55 obtained with sodium lidar 
observations is indicated with a dashed line. 

Figure 14. Horizontal wavelength versus observed wave period. The horizontal 
wavelengths were calculated using the gravity wave dispersion relation by 
assuming that the intrinsic wave period was equal to the observed period T 0 b- 
The solid line is the ML power-law fit of the form X^ = C(T ob )P where p is the 
slope and C is the coefficient. The power-law fit X z = 0.093(T ob ) 1 - 52 
obtained with sodium lidar observations is indicated with a dashed line. 

Figure 15. Vertical phase velocity versus observed wave period. The solid line is the ML 
power-law fit of the form c z = C(T ob )P where p is the slope and C is the 
coefficient. The power-law fit c z = lO-lfT^) -0 - 54 obtained with sodium lidar 
observations is indicated with a dashed line [Gardner and Voelz, 1987]. 

Figure 16. Kinetic energy versus vertical wave number. The solid line is the ML power- 
law fit of the form KE = Qlc^P where p is the slope and C is the coefficient 
The kinetic energy dependence predicted by the linear saturation theory is 
shown in the dashed line. 
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Figure 17. Kinetic energy versus horizontal wave number. The horizontal wave number 
was calculated using the gravity wave dispersion relation by assuming that the 
intrinsic wave period was equal to the observed period T 0 jj. The solid line is 
the ML power-law fit of the form KE = C(k x )P where p is the slope and C is 
the coefficient. The power-law fit KE = 5.3 x l(H(k x )* l -0 5 obtained with 
sodium lidar observations is indicated with a dashed line [Gardner and Voelz, 
1987]. 

Figure 18. Kinetic energy versus observed wave period. The solid line is the ML power- 
law fit of the form KE = C(f 0 b)P where p is the slope and C is the coefficient. 
The power-law fit KE = 9 x 10' 5 (f O b)' 159 obtained with sodium lidar 
observations is indicated with a dashed line [Gardner and Voelz, 1987]. 

Figure 19. Altitude variations of rms wind perturbations for 8 March 1984 Urbana data. 

Each 60-minute density perturbation profile was spatially low-pass filtered 
using a cutoff wavelength of 5 km. 

Figure 20. Altitude variations of rms wind perturbations averaged over eight Urbana 
observation nights. Each 60-minute density perturbation profile was spatially 
low-pass filtered using a cutoff wavelength of 5 km. 

Figure 21. Altitude profile of kinetic energy density computed from rms wind perturbation 
data averaged over eight Urbana observation nights. 

Figure A.l. Predicted shot noise floor as a function of integration time for observations 
over 35 - 50 km altitude range using the UIUC lidar system. 
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Figure A.2. Predicted vertical wavelength resolution as a function of integration time for 
observations over 35 - 50 km altitude range using the UIUC lidar system. 
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Figure 3. Spatially filtered density perturbation profiles measured on 13 August 1984 
using the UIUC lidar at Urbana, IL. The low-pass filter cutoff wave number 
was 2 x KH corresponding to a cutoff wavelength of 5 km. The diagonal 
lines indicate the apparent 0.32 m/s downward phase progression of the 7.7 
km wave. 



47 



LOCAL STANDARD TIME 


Figure 4. Spatially filtered density perturbation profiles measured on 29 February - 1 
March 1984 using the UIUC lidar at Urbana, IL. The low-pass filter cutoff 
wave number was 5.5 x 10*^ corresponding to a cutoff wavelength of 1.8 
km. The diagonal lines indicate the apparent 0.25 m/s downward phase 
progression of the 2.4 km wave. 
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Figure 5. Average vertical wave number spectrum measured at Urbana, IL on 20 
February, 8 March, 17 July, 13 August, 30 September, and 12 November 
1984 using the UIUC lidar. The altitude range of observations is 35 - 50 km 
with a spatial resolution of 300 m and a temporal resolution of 60 minutes. 
The dashed line is the estimated shot noise floor. 
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Figure 9. Seasonal distribution of horizontal wavelengths. The horizontal wavelengths 
were calculated using the gravity wave dispersion relation by assuming that 
the intrinsic wave period was equal to the observed period T ob . 
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Figure 12. Vertical phase velocity versus vertical wavelength. The solid line is the ML 
power-law fit of the form c z = C(X Z )P where p is the slope and C is the 
coefficient. The power-law fit c z = 4.09(X Z )* 0 -" obtained with sodium lidar 
observations is indicated with a dashed line [Gardner and Voelz, 1987]. 





T ob (min) 

Figure 13. Vertical wavelength versus observed wave period. The solid line is the ML 
power-law fit of the form X z - C(T ob )P where p is the slope and C is the 
coefficient, the power-law fit \ = 0.4(T ob )-°- 55 obtained with sodium lidar 
observations is indicated with a dashed line. 
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Figure 14. Horizontal wavelength versus observed wave period. The horizontal 
wavelengths were calculated using the gravity wave dispersion relation by 
assuming that the intrinsic wave period was equal to the observed period T ob . 
The solid line is the ML power-law fit of the form Xj, = C(T ob )P where p is the 
slope and C is the coefficient. The power-law fit X z = O.OQSCTqJj) 1 - 52 
obtained with sodium lidar observations is indicated with a dashed line. 
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Figure 15. Vertical phase velocity versus observed wave period. The solid line is the ML 
power-law fit of the form c z = C(T ob )P where p is the slope and C is the 
coefficient. The power-law fit c 2 = lO-lCT^)- 0 - 54 obtained with sodium lidar 
observations is indicated with a dashed line [Gardner and Voelz, 1987]. 
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Figure 16. Kinetic energy versus vertical wave number. The solid line is the ML power- 
law fit of the form KE = QkJP where p is the slope and C is the coefficient 

The kinetic energy dependence predicted by the linear saturation theory is 
shown in the dashed line. 




KINETIC ENERGY (m 2 /s 2 


60 



Figure 17. Kinetic energy versus horizontal wave number. The horizontal wave number 
was calculated using the gravity wave dispersion relation by assuming that the 
intrinsic wave period was equal to the observed period T ob . The solid line is 
the ML power-law fit of the form KE = C(k x )P where p is the slope and C is 
the coefficient The power-law fit KE = 5.3 x l<H(k x )-L05 obtained with 

sodium lidar observations is indicated with a dashed line [Gardner and Voelz, 
1987]. 
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Figure 18. Kinetic energy versus observed wave period. The solid line is the ML power- 
law fit of the form KE = C(f 0 b)P where p is the slope and C is the coefficient. 
The power-law fit KE = 9 x 10'5(f o jj)-L59 obtained with sodium lidar 
observations is indicated with a dashed line [Gardner and Voelz, 1987]. 
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Figure 19. Altitude variations of rms wind perturbations for 8 March 1984 Urbana data. 

Each 60-minute density perturbation profile was spatially low-pass filtered 
using a cutoff wavelength of 5 km. 
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Figure 20. Altitude variations of rms wind perturbations averaged over eight Urbana 
observation nights. Each 60-minute density perturbation profile was spatially 
low-pass filtered using a cutoff wavelength of 5 km. 
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Figure 21. Altitude profile of kinetic energy density computed from rms wind perturbation 
data averaged over eight Urbana observation nights. 
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Figure A.2. Predicted vertical wavelength resolution as a function of integration time for 
observations over 35 - 50 km altitude range using the UIUC lidar system. 




